Study of the magnetic state of Ko.sFei 6 Se 2 
using the five-orbital Hubbard model in the Hartree-Fock approximation 
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Motivated by the recent discovery of Fe-based superconductors close to an antiferromagnetic in- 
sulator in the experimental phase diagram, here the five-orbital Hubbard model (without lattice 
distortions) is studied using the real-space Hartree-Fock approximation, employing a 10x10 Fe clus- 
ter with Fe vacancies in a yE x v5 pattern. Varying the Hubbard and Hund couplings, and at 
electronic density n=6.0, the phase diagram contains an insulating state with the same spin pat- 
tern as observed experimentally, involving 2x2 ferromagnetic plaquettes coupled with one another 
antiferromagnetically. The presence of local FM tendencies is in qualitative agreement with Lanc- 
zos results for the three-orbital model also reported here. The magnetic moment ~3/is/Fe is in 
good agreement with experiments. Several other phases are also stabilized in the phase diagram, in 
agreement with recent calculations using phenomenological models. 



Introduction. Among the most recent exciting de- 
velopments in the field of Fe-based superconductors 1 - is 
the discovery of superconductivity (SC) with T c ~30 K 
in the heavily electron-doped 122 iron-chalcogcnidcs 
Ko.8Fe2_2;Se2 and (Tl,K)Fe2_ x Se2 compounds^ These 
materials contain ordered Fe vacancies in the FeSe lay- 
ers, increasing the complexity of these systems. Re- 
cent neutron scattering results for the parent compound 
Ko.sFei.eSean with the Fe-vacancies arranged in a \/5 x 
a/5 pattern, revealed an unexpected magnetic and insu- 
lating state involving 2x2 Fe plaquettes that have their 
four Fe spins ferromagnetically ordered, and with these 
plaquettes coupled to each other antiferromagnetically^ 
The ordered magnetic moment is 3.31 //.#/Fe. the largest 
among all Fe pnictidc and chalcogcnide superconductors, 
and the magnetic transition occurs at a high temper- 
ature T/v~559 K. Angle-resolved photoemission experi- 
ments for (Tl, K)Fei.78Se2 have revealed a Fermi sur- 
face (FS) with only electron-like pockets at the (it, 0) 
and (0, 7r) points and a nodeless superconducting gap 
at those pockets^ The superconducting phase in these 
compounds cannot be explained by the nesting between 
hole and electron pockets££ Moreover, the resistivity of 
these materials displays a behavior corresponding to an 
insulator in a robust range of the Fe concentration x£ 
suggesting that SC may arise from the doping of a Mott 
insulator, as in the cuprates. All these results certainly 
have challenged prevailing ideas for the origin of SC in 
these materials that were originally based on a nested FS 
picture and a metallic parent state. 

Several theoretical efforts have recently addressed the 
exotic magnetic state that appears in the presence of 
vacancies. Band structure calculations described this 
state as an antiferromagnetic insulator with a gap ~0.4- 
0.6 eV.— £ Several model Hamiltonian calculations have 
also been presented and, in particular, two recent pub- 
lications are important to compare our results against. 



Yu et alJ£ analyzed this problem using a phenomeno- 
logical Ji- J2 spin model (see also Ref. @) with nearest- 
neighbors (NN) and next-NN terms superexchange cou- 
plings, studied via classical Monte Carlo. In this anal- 
ysis the couplings inside the 2x2 plaquettes and those 
between plaquettes were allowed to be different, and also 
to take positive or negative values. Five antiferromag- 
netic phases, including the phase found experimentally^ 
in K0.8Fe1.6Sc2, which was dubbed "AF1", were found 
varying the J\ and J2 couplings^ From a different per- 
spective that relies on a two-orbital (d xz and d yz ) spin- 
fermion model for pnictides, and with tetramer lattice 
distortion incorporated, Yin et al^ studied the regime 
of electronic density n=l (one electron per Fe), where 
they also reported the presence of an AF1 state, found 
competing with a "C" type state with wavector (w, 0). 

In the present publication, a more fundamental five- 
orbital Hubbard model, without lattice distortions, is in- 
vestigated. Our main result is that increasing the Hub- 
bard coupling U and the Hund coupling J, a robust re- 
gion of stability of the AF1 state is found. Our effort 
allows to display the regions of dominance of the many 
competing states in terms of U and J/U, facilitating 
a discussion on possible phase transitions among these 
states by varying experimental parameters. A sketch of 
the AF1 state and its two main competitors, the C and 
AF4 states, is in Fig. [TJ Our results agree qualitatively 
in several respects with the phenomenological studies of 
Refs. [I0IIT2I particularly if a combination of results of 
these investigations is made. Finally, also note that a re- 
cent studyi^ of the three-orbital Hubbard modeli 4 - using 
mean-field techniques^ has also reported the existence of 
an AF1 state but with orbital order (OO). The relation 
with our results will also be discussed below. 

Models and methods. In this manuscript, the standard 
multiorbital Hubbard model will be used. This model 
has been extensively described in several previous pub- 
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FIG. 1: (Color online) (a) Sketch of the AF1 state found to be 
stable in a region of the U-J/U phase diagram (see Fig. 2) in 
our HF approximation to the five-orbital Hubbard model, in 
agreement with neutron diffraction.— (b) A competing state 
dubbed AF4 (stable at smaller J/U's in Fig. 2). (c) The 
C competing state. For (b) and (c), a subset of the 10x10 
cluster used is shown. 



lications, by our group and others. More specifically, 
the model used is the five-orbital Hubbard model defined 
explicitly in Ref. [TH with the hopping amplitudes intro- 
duced by Graser et alr^- By construction, this model has 
a FS that is in close agreement with band structure calcu- 
lations and angle-resolved photoemission results for the 
pnictides without vacancies. The presence of the real- 
istic AF1 state in our results, as shown below, suggests 
that the same set of hopping amplitudes can be used in 
a system with Fe vacancies. The electronic density will 
be n=6.0, i.e. 6 electrons per Fe, for all the five-orbital 
model results presented below. The couplings are the 
on-site Hubbard repulsion U at the same orbital and the 
on-site Hund coupling J. The on-site inter-orbital repul- 
sion U' satisfies U'=U-2J. The computational method 
that is employed to extract information from this five- 
orbital model relies on the study of a 10x10 cluster, as 
sketched in Fig.QJa), using periodic boundary conditions. 
In this cluster, several vacancies and 2x2 building blocks 
fit comfortably inside, giving us confidence that the main 
local tendencies to magnetic order are not dramatically 
affected by size effects. 

With regards to the actual many-body technique used 
to study the 10x10 cluster, here the real-space Hartree- 
Fock (HF) approximation was employed. The method is 
a straightforward gen eralization of that used recently by 
our group in Ref. [l7| in the study of charge stripe ten- 
dencies for the two-orbital model. This HF real-space 
approach was preferred over a momentum-space proce- 
dure in order to allow for the system to select spon- 
taneously the state that minimizes the HF energy, at 
least for the finite cluster here employed. In practice, 
the many fcrmionic expectation values that appear in 
the HF formalism must be found iteratively by energy 



minimization. At the beginning of the iterative process, 
both random initial conditions as well as initial ordered 
states favoring the many phases that are anticipated to 
be in competition were employed. After each of the com- 
puter runs using different initial conditions have reached 
convergence, at a fixed U and J/U, a mere comparison 
of energies allowed us to find the ground state for those 
particular couplings. In our setup, typical running times 
for one set of couplings U-J/U required approximately 
20 hours of CPU time to reach convergence^ Dozens of 
computer cluster nodes have been used to complete our 
analysis in a parallel manner. 
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FIG. 2: (Color online) Phase diagram of the five-orbital Hub- 
bard model with y/Exy/E Fe vacancies studied via the real- 
space HF approximation to a 10x10 cluster, employing the 
procedure for convergence described in the text. With in- 
creasing U, clear tendencies toward magnetic states are de- 
veloped. The realistic AF1 state found in neutron scattering 
experiments^ appears here above J/f7=0.15 and for U larger 
than 2.5 eV. The notation for the most important states is ex- 
plained in Fig. [T] and for the rest in Refs. Ia ficilfl2l The region 
with low-intensity yellow circles at small U is non-magnetic— 

Results. The main results arising from the computa- 
tional minimization process just described are summa- 
rized in the phase diagram shown in Fig. O Since the 
hopping parameters of Ref. [l6j are already in eV units, 
our Hubbard coupling U is also displayed in the same 
units. The notation for the many competing phases used 
here is that of Refs. to facilitate comparisons. 

The main result of the present work is that our phase di- 
agram displays a robust region where the magnetic order 
unveiled by neutron diffraction^ see Fig. QJa) , is found 
to be stable. The ratio J/U needed for the AF1 phase 
to be the ground state is in good agreement with previ- 
ous estimations for the same model, although obtained 
in the absence of vacancies, based on the comparison of 
Hubbard model results against neutron and photoemis- 
sion dataJ^ The ratio J/U is surprisingly similar between 
the pnictides and the chalcogcnides. With regards to the 
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actual value of U in eV's, the range unveiled in previous 
investigations that focused on the "1111" and "122" fam- 
ilies of pnictides was approximately 1.5 eV (see Fig. 13 
of Ref. Il5f ). The increase to 2.5 eV in the present in- 
vestigation is not surprising in view of the more insu- 
lating characteristics of materials such as K0.sFe1.6Se2, 
and suggests that merely adding vacancies to the inter- 
mediate U state of the pnictides (without vacancies) is 
not sufficient to stabilize the AF1 state but an increase 
in U is also needed. Finally, with regards to 00, none 
is observed in the AF1 state in the range of U shown 
in Fig. [5J i.e. for U<6 eV. In this range, the electronic 
density of all the orbitals (d xz and d yz in particular) is 
independent of the site location in the cluster analyzed. 
However, upon further increasing U to 8 eV and beyond, 
the same 00 pattern found in the three-orbital models 
appears in our calculations (not shown explicitly), with 
the populations of the d xz and d yz orbitals now being 
different at all sites. It seems that with five orbitals the 
AF1 state manifests itself both with and without 00, 
depending on U, while for three orbitals the intermedi- 
ate phase with AF1 magnetic order and without 00 is 
not present^ 
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FIG. 3: (Color online) Density of states of the AF1 and 
C phases sketched in Figs. [T{a,c), a t the U's indicated, 
J/{7=0.25, and using a 10x10 cluster. The gap at the chem- 
ical potential suggests that the AF1 state ({7=3 and 5) is an 
insulator, although with a mild U dependence in the value of 
this gap. On the other hand, the C state appears to have only 
a pseudogap at the Fermi level.— 

Together with the realistic AF1 phase, Fig. [2] reveals 
several other states, and two of them are prominent. 
Keeping the ratio J/U constant but reducing U, the pre- 
viously described C-type state (Fig. [TJc)) was found to 
be stable. This is reasonable since without Fe vacancies 
this state is the dominant spin order in the intermediate 
range of couplings, where the ground state is both metal- 
lic and magnetic^ In K0.sFe1.6Se2, as the bandwidth is 
increased by, e.g., increasing the pressure, a transition 
from the AF1 to the C-statc could be experimentally 
observed. In these regards, our conclusions agree with 
Ref. [H that the C-state is the main competitor of the 



AF1 state. However, note that other states reported in 
Ref. [13 are also present in our phase diagram. For in- 
stance, the AF4 state (Fig. [21(b)) is stable in a large re- 
gion of parameter space at small values of J/U. Thus, 
overall our results suppo rt a combination of the main 
conclusions of Refs. llOUl2l . 

The density-of-state (DOS) for the AF1 phase is shown 
in Fig. [3] for representative couplings. The presence of 
a gap at the chemical potential indicates an insulating 
state, in agreement with experiments.— This is not sur- 
prising considering that the transport of charge from each 
2x2 building block to a NN block may be suppressed 
due to the effective antifcrromagnctic coupling between 
blocks, at least at large U and J. In other words, us- 
ing a tilted square lattice made out of 2x2 superspin 
blocks, the state is actually a staggered antiferromagnet 
that is known to have low conductance. On the other 
hand, it is interesting to observe that the AF1 gap is 
only weakly dependent on U, suggesting that not only 
the increase in U is responsible for the insulating be- 
havior but there must be other geometrical reasons that 
may contribute to the gap through quantum interference. 
This is reminiscent of results reported years ago for the 
insulating CE phase of half-doped manganites, state that 
is stabilized in the phase diagram even in the absence of 
electron-phonon coupling due to the peculiar geometry 
of the zigzag chains involved in the CE state and the 
multi-orbital nature of the problem, that induces a band 
insulating behavior.— Thus, in agreement with recent in- 
dependent observations,— our results suggest that the 
insulator stabilized in the presence of Fe vacancies may 
have a dual Mott and band-insulating character. Note 
also that the competing C-state only has a pseudogap 
(Fig. [3]), and thus it may be a bad metal.— 

With regards to the strength of the FM tendencies 
in each of the 2x2 building blocks of the AF1 state, 
examples of the values of the magnetic moment m (in 
Bohr magnetons, assuming g=2, and at J/L7=0.25) 
are m=3.87 (17=3.0), m=3.93 (E/=4.0), and m=3.95 
({7=5.0), in good agreement with neutron diffraction 
results^ 77i=3.3. Thus, the Fe spins in the AF1 su- 
perblocks are near the saturation value 4.0 [Ib at n=6.0. 
Note that the competing C-phase also has a surprisingly 
large moment to=3.5 at {7=2.0 and J/{7=0.25. 

Results for the three-orbital Hubbard model. The re- 
sults reported thus far have been obtained under the HF 
approximation. Better unbiased approximations for this 
model are not currently available. However, at least con- 
sistency checks of the present results can be carried out 
using the Lanczos technique restricted to the 2x2 clus- 
ter of irons that forms the AF1 state. For our problem, 
an additional simplification from five to three orbitals 
(d X z j d yz , and d xy ) is needed to reduce the Hilbert space 
to a reasonable size, thus here the model introduced by 
Daghofer et al^ was used. The present Lanczos study 
is equivalent to a 12-sites one-orbital Hubbard model 
which can be done comfortably with present day com- 
puters even with the open boundary conditions (OBC) 
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FIG. 4: (Color online) Total spin of a 2x2 cluster using the 
three-orbital Hubbard model^ and the Lanczos technique, 
varying U and J/U. Panel (a) is for 10 electrons, while (b) 
is for 14 electrons. The "undoped" limit is for 16 electrons in 
the three-orbital model.— The figure shows that in these hole 
doped clusters FM tendencies develop as U grows, at realistic 
J/U's, compatible with the five-orbital HF results. 

employed tiered Our focus has been on the total spin 
quantum number to search for indications of FM tenden- 
cies in the 2x2 cluster. The main results are in Fig. 2] for 
the case of hole doping. These results indicate that with 
increasing U and J/U, FM tendencies indeed develop, 
in agreement with the five-orbital HF results (Fig. [2]). 
With electron doping also FM tendencies were found (not 
shown) . In these Lanczos results the transition from low- 
to high-spin is rather abrupt. However, note that in the 
2x2 cluster there is obviously no room to distinguish a 



fully FM state from an AF1 state, thus the large-spin re- 
gion of the 2x2 cluster may correspond to any of the two 
if larger clusters could be studied. Nevertheless, solving 
exactly this case allows us to confirm that with increas- 
ing U (at realistic J//7's), there is a clear tendency in the 
multiorbital Hubbard model towards local FM order in 
2x2 clusters. 

Summary. Real-space HF- approximation results for 
the five-orbital Hubbard model, supplemented by Lanc- 
zos calculations for three-orbital on a 2 x 2 plaquette, have 
been presented for the case of a \/T) x \/5 arrangement of 
Fe vacancies. The phase diagram obtained by varying U 
and J/U contains the magnetic state found in neutron 
diffraction experiments^ This state arises at intermedi- 
ate couplings U and J/U, and in the phase diagram it 
is not in contact with the paramagnetic metallic state of 
the weak coupling limit. Thus, FS nesting cannot explain 
the stability of the AF1 magnetic state in the presence of 
Fe vacancies^ The density of states shows that the AF1 
state is an insulator, but since the gap does not present 
a strong dependence on U its origin may reside in a 
combination of Hubbard and band-insulator features. In 
agreement with recent spir*i£ and spin-fermionaS model 
calculations, several other magnetic phases were found 
here, suggesting that transitions among these competing 
states, or among AF1 with and without OO^ could be 
observed experimentally particularly by modifications in 
the carrier's bandwidth. 
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